Conservation law of operator current in open quantum systems 
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We derive a fundamental conservation law of operator current for master equations describing 
reduced quantum systems. If this law is broken, the temporal integral of the current operator of an 
arbitrary system observable does not yield in general the change of that observable in the evolution. 
We study Lindblad-type master equations as examples and prove that the application of the secular 
approximation during their derivation results in a violation of the conservation law. We show that 
generally any violation of the law leads to artificial corrections to the complete quantum dynamics, 
thus questioning the accuracy of the particular master equation. 
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I. INTRODUCTION 

Quantum master equations are a valuable tool when 
describing the dynamics of open systems. However, the 
typically employed reduced-density-operator theory does 
not a priori guarantee that the resulting evolution main- 
tains all necessary physical properties. A well-known ex- 
ample of the pursuit for these properties is given in the 
case of quantum Markov processes by the Lindblad form 
describing the most general generators of the quantum 
dynamical semigroup P, Q ■ Even though this form and 
its time-dependent generalizations ensure certain criti- 
cal properties of quantum evolution they do not 
account for time-local conservation of observables. Ad- 
ditionally, many microscopic derivations of master equa- 
tions exploit the secular approximation [3| that has been 
shown to lead to non-physical behavior including non- 
conservation of electric charge • Related to the con- 
servation, the continuity equation for current has been 
studied in the coarse-grained description of the reduced 
quantum dynamics in Refs. 

In this paper, we introduce a general framework for the 
conservation law that all master equations for reduced 
quantum systems should ideally follow. The evolution 
obeying the law ensures that the temporal integral of the 
current operator of an arbitrary system observable, as 
obtained from the commutator with the Hamiltonian of 
the complete system, yields the change of that observable 
in time. In other words, the current flowing into the sys- 
tem equals the current obtained by it. As examples, we 
apply the conservation law to a few typical derivations 
of master equations leading to the Lindblad form and 
show that the secular approximation leads to nonconser- 
vation. Hence, Lindblad-type master equations do not 
intrinsically guarantee conservation for all observables. 



II. CONSERVATION OF OPERATOR 
CURRENT 

Let us consider a quantum system described by a den- 
sity operator p. We differentiate a subsystem S described 



by a reduced density operator ^5 = Tte{p}, where the 
trace is over the remaining environmental degrees of free- 
dom, and we denote a general S'-observable as G. We 
refer to the time derivate of the expectation value of the 
observable as operator current and write it as 



dt 



Tr' 



^~dt 



(1) 



The von Neumann equation -^p = —j-\H, p], results in 
the Ehrenfest theorem stating that [l^l 



d 
di 



(G) 



TT{p[G,H]} + Tr{p 



dG 
'dt 



(2) 



where H is the Hamiltonian of the total system. In order 
to relate this to the evolution of the subsystem of inter- 
est, we write the total Hamiltonian in the general form 
H = Hs ® Ie + Is ® He + Hi , where we have separated 
Hamiltonians for the system, the environment, and the 
interaction between them, respectively. Using the full 
form of the Hamiltonian results in 



I (G) = [Trs{ps[G. Hs]} + Tr{p[G, Hj]] 



(3) 



yielding our first definition for the operator current. We 
have denoted the trace over the subsystem degrees of free- 
dom by Trg. The current is comprised of three separate 
contributions. The first and third terms relate to the evo- 
lution of the closed system, and they are affected by the 
environment only through ps. The second term describes 
current induced by the interaction with the environment 
and vanishes for closed systems. 

To illustrate how decoupling of the eigenstate popu- 
lations and the coherence between them leads to non- 
physical behavior, we provide a simple example. Con- 
sider a two-level system whose Hilbert space is 7^5 = 
spandig) , |e)}), where Hs\i) — Ei\i) and inner prod- 
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ucts for an arbitrary system operator Os are defined 
as {s\Os\p) = Ofp, where s,p £ {ff, e}. Assume that 
the system starts from a fully excited state pf^ = 1 
and = = 0, and Hi ^ G ® E, where E is 



any nontrivial environment operator and G is time in- 
dependent. We assume that G is not diagonal in the 
eigenspace of the system Hamiltonian so that the sys- 
tem has a nonzero relaxation rate to the ground state. 
We consider a zero-temperature environment, and hence 
the system relaxes to the ground state and we have a 
stationary state p^g = 1 and pf^ = p^g = 0. The expec- 
tation value of an observable assumes the general form 
(G) = (Ggg-Gge)pf<, + 2Ke(pfgGgg) + Gee so that, in the 
long-time limit, the temporal change in the expectation 
value becomes A (G) = Ggg — Gee, which is nonzero for 
an almost arbitrary operator G. Equation ([3|) yields a 
current operator for G as 7^ = —ji[G, Hs] corresponding 
to the usual definition for subsystem current operators 
[l3| . Hence, we have {la) = —2bJoiim{pg^,G eg) , where 
woi = {Ee — Eg)/h, so that the integrated current be- 
comes J {Ig) dt = — 2a;oi3m(Geg/ Pgedt). Up to this 
point, the example has been on a very general level and 
no approximations on the dynamics have been invoked. 
However, if the populations and coherences decouple in 
the description of the dynamics for ps, our assumption 
of the initial state implies that p^^ = at all times. Thus 

{Iq) = at all times yielding A (G) ^ J {Ig) dt for an 
almost arbitrary G. Hence, the local conservation of the 
operator current breaks down in the sense that the cur- 
rent cannot accurately describe the temporal change of 
the observable. In the following, we formulate a general 
condition ensuring this conservation. 



Let the temporal evolution of the reduced system be 
described by a master equation as 



Jt^' 



--[Hs,ps] + D, 



(4) 



where we have separated the part relating to unitary 
evolution from the generator and D = D{ps,t) repre- 
sents a generalized dissipator; that is, it also accounts 
for any unitary contribution stemming from the system- 
environment interaction. Combining Eqs. ([T]) and ^ 
results in our second definition for the operator current: 



^ (G) = -'-Trs{ps[G, Hs]} + Trs{DG} 
+ Trs < PS— > ■ 



(5) 



dt 



Thus, we have two fundamental definitions provided by 
Eqs. ([2]) and ^ leading to a necessary and sufficient 
condition for the conservation of the operator current: 



:Tr{p[G,Hi]} = Trs{DG}. 



(6) 



This condition states that the dissipative current ob- 
tained from the master equation must be equal to the 
dissipative current related to the interaction Hamilto- 
nian and ensures the conservation of operator current 
which we define as A (G) = / {Ig) dt, where A (G) is 
the temporal change given by the master equation and 
Eq. ([3]) defines {Ig) = ^ (G). In practice, the complete 
dynamics of the total density operator can be unknown, 
and hence it is convenient to cast Eq. (6) into the form 
TrsiDG} = for all [G,i?/] = 0. This condition empha- 
sizes the fact that physical quantities which are conserved 
by the interaction Hamiltonian, have to be conserved in 
the reduced dynamics; that is, the dissipative current 
must vanish in this case. 

Note that the preceding derivation required that the 
master equations describe the system dynamics exactly. 
However, a typical derivation of a quantum master equa- 
tion involves a set of approximations resulting in an 
approximate description of the dynamics. The con- 
servation law is a valuable tool also in this case: We 
can take any master equation determining the reduced- 
system evolution and define a set of corresponding to- 
tal quantum states {|^)} as the ones satisfying ps — 
Tte{\^){'^\} [Hi- For each of these states, the evolution 
is unitary, and hence we can define an operator corre- 
sponding to the total Hamiltonian Ha = 11 + Hs ■ Thus 
the real approximate evolution, ps, corresponds to an 
exact evolution of a different system. As a consequence, 
Eq. ® yields -jrTT{\^){^\[G,Hi + Hs]} = Tts{DG}. If 
the condition in Eq. ([5]) is not obeyed naturally by the ap- 
proximate master equation, we obtain Hs ^ 0. Hence, an 
artificial effective Hamiltonian emerges in the complete 
description of the dynamics. Thus the conservation law 
provides an indicator of the reliability and accuracy of 
different approximations leading to reduced-system dy- 
namics even if the complete quantum dynamics cannot 
be solved. 

Let us return to the two-level example and apply the 
conservation law. We have [G, Hi\ = implying that 
the dissipative current vanishes. A general master equa- 
tion yields Trs{^G} = {Ggg - Gee)Dgg + 2^e{DgeGeg}, 
where we used Dgg = —Dee and Dge = -D*g stemming 
from the properties of the density operator through the 
master equation. If populations p^g and coherences 
decouple as in typical master equation approaches, the 
conservation law only holds for constant populations, 
which is a contradiction. Hence, the accuracy of the ap- 
proach is compromised as discussed above. 



III. PROPERTIES OF DISSIPATIVE CURRENT 

Let us define the most general form for the interaction 
Hamiltonian as Hj = Aa ® Ba where = A^, acts 
on the system degrees of freedom and = _B^, on the 
environment degrees of freedom. The dissipative current 



3 



on the left-hand side of Eq. ^ becomes 



(7) 



in derivations exploiting the secular approximation as it 
leads to the decoupling of populations and coherences, a 
feature which was shown above to result in nonphysical 
behavior in general. 



allowing us to reduce the conservation law to a compar- 
ison of traces over S. Formulating operators TrE{Bap} 
requires knowledge of the total system evolution and, 
hence, must be done for each system separately. How- 
ever, an adequate condition for the disappearance of the 
dissipative current, not dependent on the time evolution, 
is evident: if [G, Hj] = J2ai^^ ® ~ dissipa- 
tive current vanishes. The interaction Hamiltonian can 
always be expressed such that {Ba} forms an orthogo- 
nal basis of the environmental operator space, and hence 
the tensor product form implies that this condition is 
equivalent to [G, Aa] ~ for each system operator in the 
decomposition. 

A large range of microscopic derivations of master 
equations relies on the Born approximation stating that 
the environment is only weakly coupled to the sys- 
tem. Thus, the density matrix of the environment is 
assumed to be negligibly affected by the interaction so 
that p(t) w ps{t) ® pe- This results in TTE{Bap} ~ 
TrE{BaPE}psit) = {Ba)EPs{i) using ( )^ for the envi- 
ronment average. A noise source for which the environ- 
ment average of the perturbation vanishes for each a, an 
assumption used in a variety of derivations, leads appar- 
ently to a vanishing dissipative current on the right-hand 
side of Eq. ([6]). This would naively imply that any deriva- 
tion of the quantum master equation utilizing the Born 
approximation and the preceeding assumption should re- 
sult in Tys{DG} — 0. However, we will show that this 
does not generally apply and that the level, at which the 
approximation is performed, is the key. Performing it in 
the derivation of the master equation as usual allows for 
weak dissipative current, whereas performing it on the 
level of Eq. ([3]) results in the artifact of total decoupling 
of the dissipative contribution. 



IV. LINDBLAD FORM AND SECULAR 
APPROXIMATION 

To connect our general theory described above to a 
few important examples, we turn our attention to quan- 
tum Markov processes [J| and study different microscopic 
derivations leading to master equations of the Lindblad 
form. The Lindblad form describes the most general 
form that the generator of a quantum dynamical semi- 
group can take, hence guaranteeing both the semigroup 
property and the properties of the dynamical map [l|, |2| ■ 
However, the form itself is an abstract construction and 
does not imply operator current conservation. Hence, mi- 
croscopic derivations leading to specific dissipators must 
be individually studied to see if they are in accordance 
with the conservation law. We are especially interested 



A. Singular-coupling limit 

Let us begin with the so-called singular-coupling limit, 
in which the coupling between the system and the envi- 
ronment is strong compared with the system Hamiltonian 
but weak compared with the bath Hamiltonian 0, IT^ . 
The master equation reads in the Schrodinger picture [3| 



- j^[Hs + Hls,Ps] 
+ Y.^{[Ap,PsAa] 

al3 



[AfiPs,Aa]), 



(8) 



which we have left in the so-called first standard form 
that can be explicitly transformed to the Lindblad form 
by a diagonalization of the Hermitian rate matrix {70^3}. 

The Lamb shift Hamiltonian is H^s = X]q/3 "^"/s^ct^/s- 
Note that 70,^3 and Sap are dependent on the Fourier 
transforms of the environment correlation functions. The 
only requirement for the correlation functions imposed by 
the derivation is sufficiently fast decay to accommodate 
the Markovian approximation. 

Let us concentrate on the special case of vanishing 
dissipative current such that [G, Hj] = 0. Since this 
implies [Aq,,G] = for each a, it suffices to study 
Hi ~ Aa <8) Ba = A B. Using the generalized 
dissipator Dgc from Eq. ([5]), we obtain 



Trs{DscG} = - ^Trs{ips[G,i]} 



'-Trs{psA[G,A]+Aps[A, G]}, 



(9) 



where S and 7 are scalar constants. Above, we utilized 
the cyclicity of the trace. This expression vanishes due 
to the commutation of A and G and, hence, the operator 
current is conserved in the case of the vanishing dissipa- 
tive current. We emphasize that even though the master 
equation was in the first standard form and utilized the 
Born-Markov approximation, the secular approximation 
was not used in its derivation. 



B. Weak-coupling limit 

Next, we study the derivation in the weak-coupling 
limit in which the secular approximation is necessary to 
achieve a Lindblad-type master equation. Again, it is 
sufficient to study interaction Hamiltonians of the form 
Hi = A <Si B. The master equation assumes in the 
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Schrodinger picture the form 

j^Ps--{[Hs+Hls,Ps] 



(10) 



C. Weak-coupling limit for adiabatically driven 
systems 



where Hls — X^^^ 'S'('^)^^(<^)^(<^)- The eigenoperators 
are defined as A{ijj) = Tle'-t^hw n(e)^n(e'), where 11 are 
projections to the respective eigenspaces of Hs and tlie 
sum is over all eigenvalues e and e' with a fixed uj. Note 
that the master equation is of the first standard form, 
and the parameters 7(w) and S{ijj) attain a dependence 
on the frequency difference w. We obtain 

Trs{b^cG} - ^ f - '-S{uj)1:ts{A{uj)ps[G,A\uj)]} 



lM.Trs{psA^{uj)[G, A{uj)] + A{uj)ps[A\uj), G]} 



(11) 



where Z)wc corresponds to the dissipator in Eq. (|10p . As- 
suming vanishing dissipative current due to commutation 
translates to [i, G] = Y.M{^)^G] = Y.M^{^). G] = 0, 
which does not necessarily result in a vanishing expres- 
sion in Eq. ([TT|) . However, if G commutes with all the 
eigenoperators individually, the operator current is con- 
served. One way to meet this special condition is to set 
[G, n(e)] = for every e implying that the observable 
G must be diagonal in the eigenbasis of Hs and hence 
cannot induce transitions. Again, this does not hold in 
general. 



Comparison with the singular-coupling limit points to 
problems with the secular approximation. In order to 
determine if this is the cause of the nonconservation, we 
go to an earlier stage in the derivation of the master 
equation in the weak-coupling limit. Without the secular 
approximation, the Redfield-type master equation yields 
a dissipator in the interaction picture for which 



TisiDirrGi} = J2 r(c.)e— *Trs{i(c.)ps 

UJ 

x^e-'*[it(^'),G,]} + c.c., 



(12) 



where G/ = e'-^^*Ge^*^^*, r(a;) is a specific Fourier 
transform of the environment correlation functions and 
c.c. denotes a complex conjugate of the preceding term. 
Here, the construction of the eigenoperators yields 
E^,e-'*[it(c.'),G/] = EAe'"''A^i^')e-'"'',Gi] = 

e'"''J2u'{AHoj'),G]e-''^'' = e'^^'[A,G]e-'"^' = 0. 
Hence, we retrieve the operator current conservation for 
the vanishing dissipative current if the secular approxi- 
mation is not performed. 



In our last example, a time-dependent external field 
is used to drive a weakly coupled system adiabatically. 
See Refs. for recent theoretical progress in 

this field. Using a superadiabatic master equation based 
on a perturbative expansion, it has been shown for two- 
level systems that the application of the secular approx- 
imation here results in nonconservation of the operator 
current The current was found to be conserved if 

the secular approximation was dropped. To account for 
the exact effect of the steering, we approach the problem 
utilizing a modified Floquet mode basis [l^ where the 
master equation in the Schrodinger picture is given by 



-[Hs + Hls,Ps] 
n 



7(0) 
2 

E 



{[Lq,psLI] + [LoPs,lI]) 



{[Lap, PsL^cb] + [Lapps. L^a]), 



(13) 



where Hls = Y.ap S{u:c.0mP)A^{a)AIl{P), Lo = 

Y,a'^{a)Ail{a), Lap = n(a)in(/3), and n(a;) = 
|02;(O) {4'x{t)[ denotes a projection operator to the a;th 
modified Floquet mode at time t. The parameters 
LOap denote the angular frequencies when the modified 
modes are used, and the real- valued functions ^iuJap) 
and S{uiap) relate to certain Foiuier transforms of the 
environment correlation function. Note that the rates 
and projection operators are time dependent as they de- 
scribe dynamics in the Floquet basis. The generator in 
Eq. (fT5|) is of the Lindblad form at each time instant and 
is obtained by applying the secular approximation. The 
derivation is carried out for Hi = A ^ 13 but we expect 
a similar result for a general decomposition. It turns out 
that Eq. does not necessarily result in vanishing dis- 
sipative current for [A, G] = and an arbitrary noise 
source. Similarly to the nondriven system, in the spe- 
cial case of [G, n(a)] = for every a, the commutation 
leads to vanishing dissipative current. The difference in 
this special condition compared with the nondriven case 
is that instead of the observable being diagonal in the 
eigenspace of the system Hamiltonian, it needs to be di- 
agonal in the Floquet basis at all times. 



To clarify the role of the secular approximation, we can 
rewrite the master equation without applying it. In the 
interaction picture, the resulting Redfield-type dissipator 
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rSnonsec „\-,rpa 
^driven,/ giveb 

±^ (14) 

where G/ denotes again the observable in the in- 
teraction picture and U{x) = \4>x{t)) {4>x{0)\ de- 
notes a propagator for the a;th mode. Note that 
A = J^fsp' n(/5')^n(/^) that in the interaction 
picture Ai = J^pp' ^° UHPI^UW- Hence 
j:pi,,e'fod''^^^'[UHf3')Au{(3),Gi] = [Ai,Gi] = since 

[A, G] = 0. Thus, the dissipative current vanishes indi- 
cating conservation. 

V. CONCLUSIONS 

We introduced a fundamental conservation law of oper- 
ator current in open quantum systems ensuring that the 



current flowing into the system equals the current ob- 
tained by it. For example, different Lindblad-type mas- 
ter equations stemming from the secular approximation 
were found not to obey the law. In the future, our anal- 
ysis provides a basic tool for exploring the regimes of 
validity of the different approximations employed in the 
reduced-density-operator theory for open quantum sys- 
tems. The conservation law is crucial especially in cases 
where the operator current is of great interest. 
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